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ABSTRACT 

An  approximate  method  is  presented  for  predicting  elastic 
collapse  of  complete  spherical  shells  subject  to  uniform  ex¬ 
ternal  pressure.  The  shell  contains  an  imperfection  in  the 
form  of  an  isolated  flat  spot  nnd  the  snap  through  behavior  of 
the  flat  spot  region  is  analyzed.  The  existence  of  higher 
modes  is  demonstrated  and  the  effect  of  various  choices  for 
the  stiffness  coefficients  at  the  edge  of  the  flat  spot  is 
investigated. 
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INTRODUCTION 


The  effect  of  small  imperfections  on  the  elastic  buckling 
strength  of  complete  spherical  shells  subject  to  uniform  ex¬ 
ternal  pressure  has  been  well  demonstrated  both  theoretically 
and  experimentally  in  recent  years.  Several  general  purpose 
computer  programs  have  the  capability  of  analyzing  a  shell  of 
revolution  with  arbitrary  meridian  and  one  of  these  programs  [1] 
has  been  used  to  calculate  buckling  pressures  for  spherical 
shells  with  isolated  flat  spots.  The  predicted  results  agreed 
well  with  the  experimental  results  reported  in  [2] .  An  anal¬ 
ysis  which  is  valid  for  spherical  shells  with  arbitrary  axi- 
symmetric  imperfections  is  presented  in  [3] ,  while  a  Rayleigh- 
Ritz  analysis  of  the  effect  of  an  isolated  flat  spot  is  given 
in  [4],  In  addition,  an  assessment  of  "imperfection  sensitivity" 
of  spherical  shells  is  provided  in  [5] . 

The  extension  of  any  of  these  approaches  to  include  asym¬ 
metric  imperfections,  such  as  that  associated  with  the  inter¬ 
action  of  neighboring  flat  spots,  would  encounter  serious 
computational  difficulties.  An  indication  of  the  complexities 
that  would  be  involved  is  provided  by  [6] ,  in  which  a  clamped 
shallow  cap  with  an  asymmetric  imperfection  is  analyzed.  An 
alternative  approach  to  the  development  of  numerical  solutions 
to  the  exact  equations  of  asymmetrical  large  deflection  shell 
theory,  would  be  a  search  for  approximate  solutions  which  might 
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O'  provide  aii  estimate  of  the  collapse  strength  albeit  at  the 

expense  of  an  accurate  representation  of  details  of  the  stress 
distribution.  The  first  step  in  this  search  is  described  in 
this  paper. 

Another  extension  of  the  problem  of  the  elastic  collapse 
of  the  imperfect  spherical  shell  would  be  the  inclusion  of 
local  plastic  effects  in  the  flat  spot  region.  The  importance 
of  local  yielding  on  the  collapse  strength  of  imperfect  spheri¬ 
cal  shells  has  been  studied  in  [7]  using  the  finite  element 
method.  However  an  approximate  analysis  which  does  not  require 
extensive  computer  time  would  be  a  significant  aid  in  practical 
dc  igr  work 

The  specific  probl'iL  to  be  analyzed  in  this  paper  is  the 
axisymmetric  elastic  collapse  of  externally  pressurized  complete 

spherical  shells  possessing  a  local  flat  spot.  The  flat  spot 

# 

region  is  modeled  as  a  shallow  cap  which  is  elastically  sup¬ 
ported  around  its  edge  by  the  "remainder"  of  the  shell.  This 
is  the  same  problem  studied  in  [4]  using  the  Rayleigh-Ritz  method. 
The  method  to  be  described  in  this  paper  is  an  extension  of  that 
used  in  [8]  for  the  analysis  of  clamped  spherical  caps  and 
studied  in  more  detail  in  [9].  It  was  shown  in  [9]  that  points 
on  the  load-deflection  curve  of  the  cap  can  be  determined  simply 
by  solving  a  quadratic  equation.  It  is  felt  that  this  method , 
because  of  its  inherent  simplicity,  would  be  suitable  for  later 
generalization  to  either  the  problem  of  multiple,  nonisolated, 
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flat  spots  or  the  elastic-plastic  problem  if  it  can  be  shown 
to  yield  a  satisfactory  solution  to  the  isolated  flat  spot 
elastic  problem. 

In  Ref.  [9]  it  was  shown  that  the  snap  through  analysis 
of  the  clamped  spherical  cap  exhibited  higher  mode  solutions. 

It  was  postulated  in  [9] ,  and  confirmed  independently  in  [10] , 
that  these  higher  mode  solutions  can  be  of  significance  in 
explaining  the  snap  through  behavior  of  the  clamped  cap, 
particularly  when  these  higher  modes  merge  with  the  lowest 
branch  of  the  load  deflection  curve.  The  existence  of  higher 
mode  solutions  for  the  imperfect  spherical  shell  problem  will 
be  investigated  in  this  paper. 

Also  to  be  studied  in  this  paper  is  the  question  of  the 
choice  of  stiffness  coefficients  at  the  edge  of  the  flat  spot 
region.  Since  the  flat  spot  region  will  be  studied  with  the 
use  of  an  approximate  analysis,  it  is  appropriate  to  also  consider 
the  use  of  approximate  stiffness  coefficients  to  represent  the 
behavior  of  the  "remainder"  cf  the  shell. 


DERIVATION  OF  GOVERNING  EQUATIONS 

Shown  in  Pig.  1  is  the  spherical  shell  of  nominal  radius  a 
and  thickness  h,  with  a  local  flat  spot  of  radius  RQ  >  a  ex¬ 
tending  over  a  i'jas-a  circle  of  radius  R.  The  shell  is  subject  to 
unife;-  a  external  pressure  q.  Diagrams  of  the  flat  spot  region 
and  the  remainder  of  the  sphere  are  shown  in  Pig.  2  which  also 
defines  the  positive  directions  for  the  force  and  moment  re¬ 
sultants,  Hr  and  M^,  and  the  horizontal  and  vertical  displace¬ 
ments,  u  and  w,  in  the  flat  spot  region. 

It  is  assumed  that  elastic  collapse  of  the  shell  is  caused 
by  snap  through  of  the  flat  spot,  and  that  the  large  deflections 
which  accompany  snap  through  are  confined  to  the  flat  spot 
region  identified  in  Pig.  2;  then  the  deformation  of  the  remainder 
can  be  characterized  by  linear  theory.  This  is  the  same  assump¬ 
tion  made  in  [4] .  Evidence  of  the  validity  of  this  assumption 
is  provided  in  11]  and  [3]  where  plots  of  the  displacement  as 
obtained  by  numerical  solution  of  the  exact  nonlinear  equations 
show  that,  even  at  the  pressures  corresponding  to  snap  through 
collapse,  the  displacements  of  the  remainder  are  within  the 
range  of  linear  theory. 

The  total  potential  energy  can  be  expressed  as 

n  *  ni  +  n2  (1) 

where  Ii^  represents  the  total  potential  energy  of  the  flat  spot, 
{j  and  It2  is  the  total  potential  energy  of  the  remainder. 
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For  symmetrical  deformations  of  the  flat  spot,  considered 
as  a  shallow  spherical  phell  subject  to  uniform  lateral  pressure  q, 
the  total  potential  energy  is 

n  *  /*  {[  e  2  -  2 (l-vje,]  (i|)  +  (V2w)“ 

1  12 (1-v  /  0  1  *  \C 

-  11*121  ^  }  rdr  -  2irq  /*  wrdr  (2) 

r  ar  dr^  0 

where 


-  _  du  ,u  .1  >dwx  ,  r  dw 
el3r  +  r  +  2  (3r)  +  3r 


(3) 

(4) 


are  respectively,  the  first  and  second  invariants  of  the 
middle  surface. strains  and 


+ 


1  dw 
r  3? 


(5) 


The  approximate  method  to  be  used  in  this  paper  involves 

setting  e2  =  0  in  equation  (2).  This  assumption  was  used  by 

Nash  and  Modeer  in  [8]  for  the  analysis  of  clamped  spherical 

caps,  and  was  interpreted  by  Nash  and  Modeer  as  being  equi~ 

valent  to  the  assumption  that  the  radial  membrane  force  in 

the  shell  is  constant.  The  snap  through  pressure  of  clamped  caps, 

computed  using  this  assumption,  was  shown,  in  [9]  to  be  within 

about  15%  of  the  exact  values  as  long  as  the  geometric  parameter 
2  */4  1/2 

X  ■  2f3(l-v  )]  R/(2RQh)  was  less  than  X  «  5.  Since 
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„ ,  ^  *£**  ’iWf '  -‘ **' '  -■  y?  V 


Bushnell  [11]  has  shown  (see  his  Figure  6)  that  for  X  >  5,  the 
imperfect  sphere  will  fail  due  to  bifurcation  of  the  flat  spot 
region  into  the  asymmetric  mode,  it  is  felt  that  setting 
e2  **  0  for  the  elastically  supported  cap  (flat  spot  region)  will 
introduce  no  more  than  the  15%  error  obtained  in  the  case  of 
the  clamped  cap.  In  fact,  the  relaxation  of  the  clamped  con¬ 
ditions  at  the  edge  of  the  cap  is  likely  to  make  the  assumption 
of  constant  radial  membrane  force  a  more  realistic  approximation 
than  it  is  for  the  clamped  cap.  Figure  5  in  Reference  [1]  lends 
support  to  this  argument. 

Denoting  the  displacements  at  the  intersection  of  the  flat 
spot  region  with  the  remainder  by  uR  and  wR  (see  Figures  1  and  2) , 
the  total  potential  energy  of  the  remainder  can  be  written  as  [4] 

h  -  t  ?  «*  <I?>R  -  I  Hb  ur  )  (6) 

If  the  following  relations  are  used 


Bk  =  \  +  \ 


UR  “  \  *  \ 


where  the  subscripts  1  refer  to  membrane  quantities  and  the 
subscripts  2  denote  the  effects  of  bending,  then  equation  (6) 
can  be  written  as 


*2  -  J¥R  1  I  "R  ‘ll’j,  ‘  I  \  \  '  \  \ 

■  i  %  \ 1 


. . .  . . _  _ _ .  — 

"-!'^+- :  K  i. ■  ■ !  ■„ !  J  AJ^W'2  -*-  -''>*  "'  ■■&■ .  :  •  ’  •  ■  -1 


U 


where  the  principle  of  reciprocity  has  been  used  to  set 


?  \  \ 


J  hr  u 


1  R2 


.✓s 


Assuming  that  linear  shell  theory  is  adequate  for  describing 

the  behavior  of  the  remainder ,  K,?  and  can  be  expressed  in 

2 

terms  of  stiffness  coefficients  as 


\  -  KU  V,  +  K12  <a?>  R 


*Rm  K21  UR2  +  K22  <i?)R 


From  membrane  theory,  we  have 


(9) 


H. 


.  .ff-OT1,.  -  qa/ 2 


u„  =  -  qa(l-v)R/Eh 
R1 


(10) 

(11) 


Utilizing  equations  (2),  (9),  (10)  and  (11)  together  with 
the  assumption  that  e2  =  0  allows  equation  (1)  to  be  written  as 


irEh" 


R  12e, 


„  _  «»»'  r  r  1  ,  fv2--'  2  (1-v)  dw  d‘w  , 

n  “  I2TT-V*)  /„  1  -  +  <v  w)  ~ r - a?  1  rdr 

-  2”9  \  Wrdr  +  **  t  K22  lW  +  (K21-  K12>  <£>r  “r2 


2.2 


-  Ku  (UR2!  2  - 


+  qa  Up  ] 


(12) 
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Setting  the  first  variation  of  equation  (12)  with  respect  to 
u  and  w  equal  to  zero  yields  the  two  governing  differential 
equations  [8] 

ex  «  -  02h2/12  (13) 

V4w  +  g2y2w  *  a  -  2|^  (14) 


A 


where  8  is  a  constant  of  integration.  The  associated  con¬ 
ditions  at  r  *  R  are 


w_ 


where 


el  +  ^ir1  1  (K21-  K12>  <^r  -  K11  “Hj  +  901 
*£f>B  +  vta?>R+  I  |k22  0s  +  i  (K2i 

-  K12>  “Hj1  -  0 


(15) 
0  (16) 


(17) 


D  «  Eh3/12 (1-V2) . 


SOLUTION 

The  general  solution  to  equation  (14)  is  [8] 

*  -  Voter)  ♦  c2  +  ^J(|  -  *j£>  (18) 

where  those  solutions  which  are  unbounded  at  r  *  0  have  been 
suppressed  and  and  C2  are  constants  of  integration.  Then 
u  is  determined  from  equations  (3) ,  (13)  and. (18)  as 
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u  * 


■>2.2 


“  lx  {cl2  JJi2ter)  -  J0Or)J2lSr)J 


CX<J 


*A2(w*fnr4-%  >  + 


3‘d 


163 


83% 


r"  -  ad, 
D  R0 


(19) 


J0#  Jl'  J2  aPPearin9  in  the  above  equations  are  the  Bessel 
functions  of  zeroth/  first  and  second  order.  It  is  con¬ 
venient  at  this  time  to  non-dimensionalize  the  equations  by 
introducing  the  following  parameters 


P  *  r/R 
w  *  w/h 


.  ✓3(l-vt)a2q 
2Eh2 


x  «  3R 

0"  *  12(l-v2)R4/a2h2 


[Note:  A* 


uR/h' 


04  (f-)  1 

R0 


K*,  * 


11  E[12(l-v2)]1/4  <h)  *'U 

*  .  i«u-»V/4  (a, 1/2  K 


1/2  ’ 


K, 


'22 


Eh 


22 


K*  *  112(1-v^)11/48  (aX/2 
*12  Eh  K 


12 


r*  «  [12(l-v2)]X/4e  .a, 1/2 


'21 


Eh 


9  K 


21 


C1/h 


Co/h 


(20) 
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■% 

.  A 

;  i 


»  i 


i  * 


Then  equations  (18)  and  (19)  transform  to 

w  *  c[  JQ(xp)  +  C*  +  - p2' *1 -  (f-)) 

10  2  2 [12 (1-v2)] 1/2  x2  R0 

2  2  2 

y  *  I  {~  12  "  C1  T~  IJi2<xP>-  J0(xp)J2(xp)J 

464Q  J2(xp)C*  e4p2  2 

x2  [12 (1-v2)  J1^  48 (1-v2)  Ro 

-  *4^1) 


The  boundary  conditions  (15)  -  (17)  become ,  in  nondimensionalized 
form 


/  n<»/i .2\  1 1/2  „*  x  dw  12 (1-v2)  , 

{[12  (1-V  )]  (K21-  K12)  g- - 9 - KX1U  }  psl 


X2  -  4Q82 


,d2w  .  ..  dw  A  av*  dw  .  [12 (1-v2)] 1/2  v*  x..  1  _  «  /ocx 

{TT  +  v  3p  +  6K22  3p  +  - 2 - (K21”  K12)ii2>  “  0  (25) 


The  imposition  of  the  first  of  these  boundary  conditions  leads 

* 

directly  to  the  following  expression  for  C2 


C2  «  -  Cx  JQ(x)  - 


2 [12 (1-v*) J 


-  (|-)] 

x2  Ro 
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Substituting  equations  (21),  (22)  and  the  nondimension- 

alized  form  of  (11)  into  the  remaining  two  boundary  conditions 

S  * 

(24)  and  (25)  results  in 


»!«:*) 2  +  (B2+B3Q)  cj  +  Djl  +  D2Q  +  D3Q2  =  o 


,  (2.7) 


B4(C*)2  +  (B5+B6Q)  c*  +  d4  +  D5Q  +  D6Q2  =0  *•.  1  (2( 

« 

where  the  B^.^Bg,  D^.^.Dg  coefficients  do  not  explicitly  . 

involve  the  pressure  parameter  Q  and  are  listed  below: 

3/2 

Bi  ■  ‘4<x>  -  J0(X>J2<X>J  *2  4 


B2  - 


2(1-v,<J)xJ1(x)K12 


*3  -  V*>  4 


B4  -  tJi2<x>  '  J0(x)J2(x))x2  K*2 


2  V4 

B5  b  •I1-1-e-V-—  Kl-vJJ^x)  -  xJc(x)  -  0J1(x)K22]x 

2  1/4  03J,(x)k!2 

B,  *  -2  [12  (1-v  j  ]  - -  — 

6  x2 


2  ^2  2  (l-u2)K* 

[12(1-0  3  ,  -X^  _  v  )K11 

\  O  “  o 


*2 - 3 — ”  l  — ~ — (|")  "  x2  1 

0*  20J  •  4 (1-v2)  *o 


+  2(g-)  K*2  ) 


12 


;  P  l«^'<7^EV:  <v  7.  JSgwSsT 


V^h-^k-.^x', 
2  1/2 

n  n  112  U-V5)]  n5  * 

D3  "  ^7T  9  K11 


24x 
-  1 


(12U-V*)) 


17* 


{ (g-)  e  (x+v)  +  ■i±2Ti-K12  [ 


(1-v2)  *  r  e4  ,a  2 


4 (1-v2)  Ro 


(£->  - 


x2] 


Dg  - 


°6  ~ 


26 


2  1/4 
[12 (1-v^) 3 


+  ea  <§-)  k*2  } 

o 

t  2fo±vI  .  (1-v)K*2  ♦  £  k*2  I 


07k? 


12 


2x* [12 (1-v^) ] 


1/4 


(29) 


*  * 

The  reciprocity  theorem  was  used  to  set  K21  «=  -  K12  in  the 

derivation  of  the  above  expressions. 

If  equation  (27)  is  multiplied  by  B4,  equation  (28)  by 

B^,  and  then  the  two  resulting  expressions  subtracted  from 

each  other,  the  following  equation  for  is  obtained: 

B^D^DgQ+DgQ2)  -  B4(D1+D2Q+D3Q2) 

C1  “  fi^B2-BiB5)  +  (B^-BjBgJQ 


(30) 


This  is  then  substituted  back  into  equation  (27)  to  obtain, 
after  some  lengthy  algebra, 

G4Q4  +  G3Q3  +  G2Q2  +  GXQ  +  6q  -  0 

■»  •» 


(31) 


where 

G4  =  (b1d6-b4d3) 2bx  +  (b1b6-b4d3)(b4b3-b1b6;b3 
+  (b4b3-Bib6)2d3 

G3  =  2  <B3D5-B4D2) (BiD6-B4D3)  B3 
+  (b1d5-b4d2) (B4B3-BlB6)  B3 

+  (BXD6-B4D3)  t2B2B3B4-B1(B2B6+B3B5)] 

•  » 

+  D2(B4B3-BlB6)2  +  2D3(B4B2-BlB5)(B4B3-BlB6) 
g2.  =  2B1(B1D4-B4D1HB1D6-B4D3)  +  (b1d5-b4d2)2b1 
+  (b1d4-b4d1) (B4B3-BlB6)B3 
+  (BlD5-B4D2)  [2B2B3B4-B1(B2B6+83Bs)] 

+  <Bid6-b4d3> (B4B2-BlB5)B2 

+  d1(b4b3-b1b6)2  +  2d2(b4b2-b1b5xb4b3-b1b6) 

+  DjIB^j-B^j)2 

Gj  -  2B1(B1D4-B4D1)(B1D5-B4D2) 

+  (b1d4-b4d1)  128^384-8^^(328^8365)] 

+  Bo(B,De-B,Do)  (8,80-8,8*) 


— - .-.-ca 


+  2D1(B4B2-B1B5)  (b4b3-e1b6)  +  d2(b4b2-b1b5) 

G0  -  (B1D4-B4D1) 2  Bx  +  (B1D4-B4D1)  (B4B2-BiB5)B2 

+  d1(b4b2-b1b5)2  (31 

Hence/  the  equation  for  determining  Q  is  a  guartic,  as  com¬ 
pared  to  a  quadratic  equation  for  Q  in  the  case  of  the  clamped 
cap.  However,  by  utilizing  equations  (29)  it  is  found  that 


B^Dg  -  B^D3  =  0 


B4B3  -  Bl®6  5  0 


This  leads  to  G4  *  G3  =  0  so  that  equation  (31)  does  in 
fact  simplify  to  a  quadratic  in  Q 


g2Q  +  GXQ  +  Gq  =  0 


with 


G2  a  *l\  +  T1T2B3  +  T22d3 


G1  =  2B1T1T3  +  T3T2B3  +  T2  D2  +  B2T1T2 


Go  -  t32b1  +  t2t3b2  +  t22dx 


where 


T1  *  B1D5  "  B4D2 


T2  *  B4B2  "  B1B5 


T«»  *  B.Di  —  B.D, 


r 


ft 


Since  the  G2,  G^,  GQ  coefficients  are  homogeneous  in  T.^,  T2»  rf3r 
any  factor  common  to  T^,  T2,  can  be  cancelled.  Substituting 
from  equations  (29)  into  equations  (36)  yields 


*  *2  2 .  *  *  *2 
T1  =  K11  0<1+v)  “  2Ki2  x  +  6  (K11K22+  2K12  ) 


T2  “  tl2-(1~^-)  iV2  x3  <V*>  'KUK22  +  2KI22> 


K 


__li  [(l-vJJ^x)  -  X  J0(x)J) 


Ta  “  -  *2  1x2k12"  <t^>  6<1+v>Kn 


2  a  *  *  *  2 

-  6  ‘SJ*  <K11K22+  2K!2  > 1 


(37) 


5/4 

where  the  common  factor  [12(l-v2)3  [J12(x)~  JQ(x) J2(x) J/246 


has  been  cancelled. 

The  load-deflection  curve  for  a  given  shell  (6,  a/RQ,v)  is 
obtained  as  follows:  For  a  given  value  of  the  strain  parameter 
x,  the  coefficients  are  evaluated  using 

equations  (29)  and  (37) .  Then  the  coefficients  G2,G^,GQ  are 
found  from  equations  (35) .  From  equation  (34) ,  two  values  of 


the  pressure  Q  are  found.  For  each  value  of  pressure,  the  con- 


ft 

stants  and  C2  are  determined  from  equations  (26)  and  (30) . 


Note  that  as  a  consequence  of  the  relations  (33) ,  equation  (30) 
reduces  to 


*  t1q+t3 

U1  T„ 


(38) 
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With  and  C2  determined,  the  displacements  to  and  y  are  given 
as  functions  of  position  p  by  equations  (21)  and  (22) . 

In  order  to  generate  the  load-deflection  curve  it  is  con¬ 
venient  to  introduce  the  average  deflection  parameter 


2  1/  2  -I 

5,  21,12(1-v2)1 -  / 


wpdp 


Though  co  does  not  include  the  deflections  of  the  "remainder" 
of  the  shell,  it  is  still  a  suitable  parameter  for  locating 
Qcr'  the  £irst  local  maximum  on  the  Q  vs.  w  curve.  Utilizing 
equations  (21)  and  (26) ,  w  can  be  expressed  as 
2  1/2  2 

w  =  -  c*  J2(x)  -  |  -  (g-)]  (4C 

Thus,  for  each  value  of  x,  the  two  values  of  Q  and  their 

associated  values  of  w  determine  two  points  on  the  load- 

deflection  curve.  The  lowest  branch  of  the  load  deflection 

curve  is  generated  by  letting  x  increase  from  zero  to  the 

2 

value  xfl  which  makes  G1  -  4G2GQ  vanish.  Those  values  of 

2 

x  >  x^  f«">r  which  G^  -  4G2Gq  is  positive  define  points  on  the 
higher  mode  loops  (see  Ref.  [9]). 

As  in  [9]  it  is  possible  to  obtain  relatively  simple  ex¬ 
pressions  for  the  "fully  snapped  through"  state  (the  configura¬ 
tion  after  snap  through  for  which  there  is  no  middle  surface 
strain)  by  considering  the  limiting  values  of  Q  and  w  as  x  *►  0. 
The  roots  of  equation  (34)  are 

°1,2  “  1-°1  ±  <G12  -  4G2G0)1/2]/2G2  141 
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Combining  equations  (29),  (37),  and  (35),  it  is  seen  that  for 
x  <<  1,  Go  is  of  the  order  x8  while  both  and  G2  are  of  the 
order  of  x6.  Thus 


lim  Q.  = 


lim  Q_ 
x->0  * 


-G1 
lim  (-—) 

x^o  ^2 


Evaluation  of  the  limiting  value  of  the  ratio  in  the  above 
equation  leads  to 


lim  Q„  =  — 


5Kue 


It  can  also  be  shown  that  for  x  «  1,  equation  (40)  reduces  to 
e£  Xl<3+v)'fe(KilK22+  2K*22)~  32K*2/ei 

32  ,KIl<1+v>  +  9<KU*22+  2K12)i  3  Q 

so  that  utilizing  the  above  expressions  for  the  limiting  values 
of  Q,  it  is  seen  that 


lim  w,  = 


*  * 


lim  to. 


-3  [Kn<3+V)+  6(KnK22+  2K12)-  32K12/6] 

"[kJ1(1+v)+  8(^*2+  2K^2)  J 


The  points  (0^,0)^)  and  (Q2,w2)  for  x  =  0  represent  respectively 
the  no  load  and  wfully  snapped  through"  states. 
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STIFFNESS  COEFFICIENTS 


The  solution  of  the  problem  cannot  be  obtained  without 

ft  fc  ft 

knowing  the  stiffness  coefficients  K^,  K12,  K22«  Essentially 
exact  expressions  for  the  stiffness  coefficients  of  exter¬ 
nally  pressurized  spherical  shells  including  the  pressure 
effect  were  derived  by  Eushnell  [12].  Bushnell's  expressions 
for  the  stiffness  coefficients  of  the  "remainder”  of  a  spheri¬ 
cal  shell  were  simplified  by  Koga  and  Hoff  [4]  in  their  flat 
spot  analysis  for  the  case  of  an  almost  complete  "remainder." 
This  corresponds  to  assuming  that  the  flat  spot  is  small  in 
extent  (R/a  <<  1) ,  an  assumption  which  has  already  been  used 
in  equation  (10)  of  this  paper.  The  stiffness  coefficients 
used  in  [9]  are; 

Ku  =  (1~q2)  {ri2  +  h2) 


*  A*  1+0  l—O 

K12  -  inr  t  <^>  <VrIoRi)  -  ^oVW  3 

o 


it  o  1/2  5 

K22  “  “  K?  {(1"Q  3  <R0  +  V > 

+  (1-v2) (1-Q2)1/2  iRx2+  I12)/e2 


+  2[l~Q(l+v)]  (iy^)1/2  (R^q-  RoIx)/0 
1-0  1/2 

-  2[l+Q(l+v)J  (RqV  I0I13/6  3 


where 


1/2  1/2 
A*  =  (1+2Q)  (^)  (R1Rq+  I1IQ)  +  (1-2Q)  (^)  (1^-  R^) 

1/2 

+  (V-1)(1-Q2).  (Rj2  +  I12)/e 

aJ  =  (i-v2)(i-q2)  (rx2  +  i12)2/e2 

1/2 

+  2[l-Q(l+v)3  (1+Q)  (R^Ig-  RqIjl)  (rx2+  I12)/e 

1/2 

-  2[l+Q(l+v)J  (1-Q)  (i±S)  <RiV  I0I1)  (R12+  Il2)/0 
+  (1+2Q)  (1-Q)  [  (RqR1)  2  +  a0i2)}/2 
+  (1-2Q)  (1+Q)  [(IqR^)2  +  (Rg^)2!^ 

+  2QR0I()R1X1  -  (1-Q2)1'2  [R0I0  (Rx2-  I-^} 

+  R^  (I02-  Rq2)]  (46) 

Rq  and  Xq  in  the  above  equations  are  the  real  and  imaginary 
parts  respectively  of  the  Hankel  function  of  the  first  kind 
of  order  zero.  Similarly  R1  and  I1  are  the  real  and  imaginary 
par  is  of  the  Hankel  function  of  order  one.  The  argument  of 
the  Hankel  function  is 

„  1/2 

z  =  0IQ  +  i  (1-  Q*))  (47! 

If  equations  (45)  are  used  in  the  expressions  for  Blf...,T3 
[equations  (29)  and  (37)3,  it  is  readily  apparent  that  the 
coefficients  G~,  G. ,  Gn  in  equation  (35)  will  involve  the  load  Q 


in  an  extremely  complicated  fashion.  However  the  very  simple 
explicit  involvement  of  Q  in  equation  (34)  can  still  be  used 
to  advantage  in  obtaining  a  solution  as  follows:  i'irst,  B^,... 
T3  are  evaluated  by  setting  Q  =  0  in  equations  (4 5) -'(47).  Then 
when  a  value  of.Q.is  determined  from  equation  (34)  this  value 
is  used  to  re-evaluate  the  stiffness  coefficients  from  equa¬ 
tions  (45) -(47).  This,  iterative  procedure  is  continued  until 
successive  values  of  Q  agree.  Then  the  solution  procedure 
continues  with  the  determination  of  u.  There  is  of  course  no 
guarantee  that  the  iteration  process  described  above  will 
converge.  Iterative  solutions  for  clamped  caps  do  encounter 
convergence  difficulties  in  the  neighborhood  of  Q ,  and  in 
fact  the  pressure  at  which  this  loss  of  convergence  occurs  is 
often  defined  as  (see  Ref.  [9]  for  a  more  detailed  discus- 
sion  of  this.)  On  the  other  hand,  the  advantage  of  extending 
the  technique  of  [9]  to  the  subject  problem  of  this  paper  is 
that  the  complete  load-deflection  curve  (including  higher  modes) 
can  be  generated  without  recourse  to  iterative  solutions  and 
their  associated  convergence  difficulties.  If  the  iteration 
process  described  above  does  lead  to  convergence  difficulties, 
the  drastic  simplifications  due  to  the  use  of  the  approximate 
solution  for  the  flat  spot  region  are  negated.  Hence  it  is 
appropriate  to  consider  using  less  exact  stiffness  coefficients 
and  in  particular,  stiffness  coefficients  which  allow  the  load- 
deflection  curve  to  be  generated  without  resorting  to  iterative 
techniques . 
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The  reason  that  the  use  of  equations  (45)  for  the  stiffness 
coefficients  requires  iteration  to  generate  the  load-deflection 
curve  is  that  the  coefficients  in  equation  (34)  depend  upon  the 
loading  parameter  Q.  If  stiffness  coefficients  which  are  inde¬ 
pendent  of  Q  are  used  instead  of  equations  (45)  ,  then  the  coef¬ 
ficients  in  equation  (34)  will  be  independent  of  Q,  and  a 
closed  form  solution  for  Q  is  immediately  obtainable.  One  such 
set  of  stiffness  coefficients  is  obtained  by  setting  Q  =  0  in 
equations  (45) -(47).  The  resulting  stiffness  coefficients  are 
the  shallow  shell  approximation  to  the  stiffness,  coefficients 
obtained  by  Baker  and  Cline  113]  in  terms  of  Thomson  functions. 
Since  the  effect  of  the  external  pressure  on  the  stiffness  of 
the  "remainder"  is  ignored  in  these  expressions ,  the  results 
will  become  inaccurate  if  the  loading  on  the  complete  sphere 
approaches  the  buckling  load  of  the  "remainder,"  i.e.,  as  Q  1. 
Thus,  the  use  of  the  Baker  and  Cline  stiffness  coefficients 
will  not  permit  the  present  analysis  to  reduce  to  that  of  a 

perfect  sphere  as  (§-)  1.  On  the  other  hand,  the  analysis 

0  a 

will  become  more  accurate  as  (=— )  departs  further  from  unity, 

R0 

i.e.,  as  the  flat  spot  becomes  "flatter."  Such  an  analysis 

serves  to  complement  the  imperfection  sensitivity  study  of 

Hutchinson  [5]  which  loses  accuracy  as  =—  departs  from  unity. 

R0 

As  indicated  above,  the  involvement  of  Q  in  equations 
(45) -(47),  not  only  explicitly  but  also  in  the  argument  of 
the  Hankel  functions,  leads  to  concern  regarding  convergence 
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of  any  iteration  process.  This  led  to  consideration  of  a  set 
of  stiffness  coefficients  which  are  independent  of  Q.  Another 
approach  to  simplifying  the  problem  while  still  retaining  the 
effect  of  Q  on  the  stiffness  coefficients  is  to  obtain  asymptotic 
values  for  the  stiffness  coefficients  for  large  values  of  0. 

By  taking  the  limit  as  0  -*■  »,  equations  (45)  reduce  to 

K11  s  ~  t2<l-Q)] 1/2  • 


K*2  =  [2(1-Q)]1/2  (48) 

These  stiffness  coefficients  are  the  counterpart  of  those  which 
were  originally  derived  by  Nachbar  [14]  for  internally  pressurized 
spherical  shells.  Nachbar' s  original  derivation  was  based  on 
the  assumption  that  the  edge  angle  of  the  "remainder"  is  close 
to  n/2,  i.e.,  that  the  "remainder"  is  very  nearly  hemispherical. 
However  no  such  restriction  is  required  to  derive  equations  (48) 
from  equations  (45) ,  a.;^  since  equations  (45)  are  valid  for  an 
almost  complete  "remainder"  (R/a  <<  1) ,  so  are  equations  (48) . 

This  raises  an  interesting  point  regarding  the  interpretation  of 
the  parameter 

04  =  12  (1-v2)  (|)2  (|)2 

ft 

In  most  previous  work  (j^-)  is  considered  as  fixed,  so  increasing 

R 

0  is  associated  with  increasing  (-)  and  hence  a  less  shallow 

cl 

p 

"remainder."  However  it  is  clear  that  (-)  can  be  kept  fixed  and 

a 
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increasing  0  can  then  be  interpreted  as  corresponding  to  a 
thinner  "remainder."  This  latter  interpretat;  on,  which  is  used 
in  this  paper,  makes  it  possible  to  utilize  equations  (48)  for 
other  than  nearly  hemispherical  "remainders." 

Since  equations  (48)  were  derived  by  taking  the  limit  as 
0  <«>,  the  validity  of  using  these  expressions  for  the  subject 

problem  where  the  range  of  interest  is  0  <  6  has  to  be  examined. 

It  was  shown  by  Cline  115]  in  his  study  of  the  effect  of  internal 
pressure  on  the  behavior  of  spherical  shells  that  the  influence 
coefficients  rapidly  approach  their  asymptotic  values  at  relatively 
small  values  of  0.  In  fact  Cline  proposed  that  the  asymptotic 
values  of  the  influence  coefficients  be  used  for  0  >_  3.  For 
externally  pressurized  spherical  shells,  Bushnell  [12]  showed 
that  when  Q  >  0.5  the  influence  coefficients  become  infinite  at 
values  of  0  which  depend  on  Q.  However,  according  to  Bushnell, 
the  stiffness  coefficients  for  externally  pressurized  spherical 
shells  are  well  behaved  and  non-zero  for  all  0.  It  is  reasonable 
then  to  assume  that  the  stiffness  coefficients  for  externally 
pressurized  spherical  shells  possess  asymptotic  behavior  as 
0  *►  *°  which  is  similar  to  that  of  influence  coefficients  for 
internally  pressurized  spherical  shells,  i.e.,  the  asymptotic 
form  given  by  equations  (48)  are  sufficiently  accurate  when 
0  >  3. 

The  use  of  equations  (48)  instead  of  equations  (45)  repre¬ 
sents  a  drastic  simplification  in  the  form  of  the  coefficients 
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in  equation  (34) .  However  these  coefficients  will  still  involve 
terms  in  which  Q  is  raised  to  a  non-integer  power  and  a  closed 
form  solution  for  Q  will  not  be  obtainable.  While  the  iteration 
procedure  required,  to  solve  the  governing  equation  is  likely 
to  be  more  stable  than  when  equations  (45)  are  used,  there  still 
cannot  be  any  guarantee  of  convergence. 

The  asymptotic  stiffness  coefficients  can  be  simplified 
even  further  by  neglecting  the  pressure  effect  (setting  Q  =  0 
in  equations  (48)).  The  expressions  then  reduce  to  the  well 
known  Geckeler  [16]  form.  However  whereas  the  Geckeler  stiff¬ 
ness  coefficients  have  previously  been  thought  to  be  restricted 
to  nearly  hemispherical  "remainders,”  the  systematic  derivation 
of  these  expressions  from  equations  (45)  shows  that  they  are 
applicable  to  any  sufficiently  thin  "remainder."  The  use  of 
the  Geckeler  stiffness  coefficients  allows  equation  (45)  to  be 
solved  directly  via  the  quadratic  formula.  The  use  of  these 
coefficients  can  be  expected  to  lead  to  inaccuracies  when  either 

0  <  3  or  ~  1. 

R0 


RESULTS 


A  computer  program  was  written  to  evaluate  the  coefficients 
G2/  ,  Gq  of  equation  (34)  and  to  solve  equation  (34)  either 

directly  when  G2,  G^,  Gq  are  independent  of  Q  or  by  iteration 
when  G2,  G^,  Gq  involve  Q.  Solutions  were  obtained  for  the 
four  sets  of  stiffness  coefficients  described  in  the  previous 
section. 

As  expected,  the  use  of  either  equations  (45)  or  (48)  did 
lead  to  convergence  difficulties,  primarily  when  attempting 
to  compute  values  of  Q  in  the  neighborhood  of  unity.  One 
reason  for  this  can  best  be  seen  from  an  examination  of  equa¬ 
tions  (48)  which  become  imaginary  when  Q  >  1.  Hence  any  itera¬ 
tive  process  which  provides  intermediate  values  for  Q  which 
are  greater  than  unity  is  doomed  to  failure.  Since  one  of  the 
main  purposes  of  undertaking  the  analysis  described  in  this 
paper  was  to  avoid  convergence  difficulties  associated  with 
iterative  solutions  of  the  problem,  no  attempt  was  made  to 
refine  the  iteration  process  described  earlier  in  this  paper. 

The  computations  revealed  the  presence  of  higher  mode 
solutions  and  a  load-deflection  curve  for  0=6,  a /Rq  =  1/1.05 
which  includes  a  higher  mode  in  the  form  of  an  isolated  loop 
is  shown  in  Figure  3.  This  result  was  obtained  using  the  exact 
stiffness  coefficients  (equations  (45) ]  and  no  convergence 
difficulties  was  encountered.  Calculations  using  equation  (48) 
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were  also  conducted  for  this  case  and  the  results  for  the  main 
branch  of  the  curve  duplicate  those  shown  in  Figure  3.  No 
higher  mode  solutions  were  sought  with  the  use  of  equations  (48) . 

A  systematic  study  of  higher  mode  solutions  for  other  values 
of  0/  a/Rg  and  other  choices  of  the  stiffness  coefficients  was 
not  undertaken,  since  from  the  results  shown  in  Figure  3  these 
solutions  appear  to  be  quite  similar  qualitatively  to  the  higher 
mode  solutions  for  clamped  caps  [9] .  The  influence  of  these 
higher  modes  on  the  snap-through  behavior  of  the  imperfect 
spherical  shell  can  only  be  explored  by  refining  the  approximate 
solution  for  the  behavior  of  the  flat  spot  region  along  the 
lines  described  in  [17] . 

A  comparison  of  the  load-deflection  curves  for  various 
choices  of  the  stiffness  coefficients  is  shown  in  Figure  4  for 
0  =  3,  a/Rg  -  1/1.15.  Notice  that  the  effect  of  neglecting  Q 
in  the  stiffness  coefficients  raises  the  value  of  Qcr*  This  is 
to  be  expected  since  the  effect  of  the  external  pressure  on  the 
“remainder"  is  to  decrease  its  stiffness.  It  is  interesting  to 
note  that  the  asymptotic  stiffness  coefficients  [equations  (48)3 
yield  results  which  are  in  excellent  agreement  with  those  ob¬ 
tained  using  the  exact  stiffness  coefficients,  even  at  this  small 
value  of  0. 

Curves  similar  to  those  of  Figure  4  were  generated  for  vari¬ 
ous  values  of  0  and  a/Rg.  The  snap  through  pressures  obtained 
from  these  curves  are  plotted  in  Figure  5  along  with  the  numerical 
results  from  [11] .  It  is  seen  that  the  approximate  solutions 


of  this  paper  give  relatively  good  results  when  0  ~  3.  The 
errors  increase  as  0  increases,  particularly  for  the  "flatter” 
imperfections.  The  results  show  that  no  significant  advantage 
in  accuracy  can  be  gained  by  using  more  accurate  stiffness 
coefficients  for  the  remainder. 

Several  errors  in  the  original  version  of  equations  (29) 
were  found  as  this  report  was  being  prepared.  They  were  cor¬ 
rected  in  the  manuscript  so  equations  (29)  are  correct  as 
presented  in  this  report.  However  these  errors  were  discovered 
too  late  to  be  corrected  in  the  computer  program.  Hence  the 
results  shown  in  Figures  3-5  are  not  correct.  It  is  expected 
that  the  qualitative  nature  of  the  results  shown  in  Figures  3 
and  4  will  not  be  affected  by  new  computations  based  on  the 
corrected  equations;  and  that  significant  improvement  in  the 
accuracy  of  the  results  shown  in  Figure  5  will  be  achieved  as 
a  consequence  of  incorporating  the  corrections  in  the  computer. 
These  new  results  will  be  incorporated  into  the  report  prior 
to  its  release  and  distribution  according  to  the  attached  list. 

Since  the  results  indicate  that  the  use  of  the  simplest 
stiffness  coefficients  (equations  (48)  with  Q  =  0)  is  appro¬ 
priate  for  this  problem,  explicit  evaluation  of  the  fully 
snapped  through  configuration  is  possible.  Substituting  equa¬ 
tions  (48)  with  Q  =  0  into  equations  (43)  and  (44)  leads  to 

°2  C  572 

“  =  572('1+v>  (5+2v+  tH) 


lim 

x*0 

lim 

x-*0 
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It  is  of  interest  to  note  that  these  values  are  independent  of 
a/Rg  and  that  the  value  for  Q2  is  almost  identical  to  that 
obtained  in  [9]  for  the  clamped  cap. 
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